NASA  Contractor  Report  198167 
ICASE  Report  No.  95-42 


ICASE 


EXPERIMENTAL  CONFIRMATION  OF  A  PDE-BASED 
APPROACH  TO  DESIGN  OF  FEEDBACK  CONTROLS 


H.  T.  Banks 
Ralph  C.  Smith 
D.  E.  Brown 
R.  J.  Silcox 
Vern  L.  Metcalf 

19950712  033 

Contract  No.  NASl-19480 
May  1995 


Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  VA  23681-0001 


Operated  by  Universities  Space  Research  Association 


DTH  QUALITY  INSPECTED  8 


EXPERIMENTAL  CONFIRMATION  OF  A  PDE-BASED 
APPROACH  TO  DESIGN  OF  FEEDBACK  CONTROLS  ' 


H.T.  Banks 

Center  for  Research  in  Scientific  Computation 
North  Carolina  State  University 
Raleigh,  NC  27695 


Ralph  C.  Smith 
Department  of  Mathematics 
Iowa  State  University 
Ames,  I A  50011 


D.E.  Brown 

Lockheed  Engineering  and  Sciences  Company 
NASA  Langley  Research  Center 
Hampton,  VA  23681 


R.J.  Silcox 
Acoustics  Division 
NASA  Langley  Research  Center 
Hampton,  VA  23681 


Vern  L.  Metcalf 

U.S.  Army  Research  Laboratory 
NASA  Langley  Research  Center 
Hampton,  VA  23681 


ABSTRACT 

Issues  regarding  the  experimental  implementation  of  PDE-based  controllers  are  discussed  in 
this  work.  While  the  motivating  application  involves  the  reduction  of  vibration  levels  for  a  cir¬ 
cular  plate  through  excitation  of  surface-mounted  piezoceramic  patches,  the  general  techniciues 
described  here  will  extend  to  a  variety  of  applications.  The  initial  step  is  the  development  of 
a  PDE  model  which  accurately  captures  the  physics  of  the  underlying  process.  This  model  is 
then  discretized  to  yield  a  vector-valued  initial  value  problem.  Optimal  control  theory  is  used  to 
determine  continuous-time  voltages  to  the  patches,  and  the  approximations  needed  to  facilitate 
discrete  time  implementation  are  addressed.  Finally,  exj^erimental  results  demonstrating  the 
control  of  both  transient  and  steady  state  vibrations  through  these  techniques  are  presented. 
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1  Introduction 


An  increasingly  popular  method  for  controlling  structural  vibrations  is  through  the  use  of  piezo¬ 
ceramic  patches  bonded  to  or  imbedded  in  the  structure.  These  patches  exhibit  the  piezoelectric 
property  that  inplane  strains  are  generated  in  response  to  an  applied  voltage.  Depending  upon 
the  geometry  of  patch  placement,  location  with  respect  to  the  structure’s  neutral  surface,  and 
the  method  of  excitation,  this  provides  a  mechanism  for  generating  both  inplane  forces  and/or 
bending  moments  in  the  underlying  structure. 

The  advantage  of  using  such  patches  as  actuators  in  many  applications  is  due  to  the  fact  that 
they  are  lightweight,  space  efficient,  relatively  inexpensive,  and  provide  a  means  of  obtaining 
structural  control  without  significantly  changing  the  passive  structural  dynamics  (they  do  not 
mass  load  the  structure  in  the  manner  of  a  shaker  or  proof  mass  actuator).  Due  to  their 
ceramic  nature,  they  can  also  be  molded  in  a  variety  of  shapes  so  as  to  fit  the  structure  under 
consideration.  Moreover,  rigid  body  torciues  and  spillover  effects  are  minimized  due  to  the  fact 
that  they  are  fully  self  contained  and  distributed  in  nature.  Finality,  they  also  exhibit  the  inverse 
piezoelectric  property  and  hence  generate  a  voltage  in  response  to  strains  in  the  material.  Hence 
a  single  jratch  or  patch  pair  can  be  used  for  either  sensing  or  actuation.  This  contributes  to 
their  efficiency  and  application  in  “smart  material”  structures. 

In  addition  to  their  utility  in  purely  structural  applications,  they  are  also  finding  increasing 
use  in  structural  acoustic  and  fluid/structure  applications.  Again,  their  advantage  lies  in  the  fact 
that  they  provide  an  efficient  means  of  controlling  the  structure  without  significantly  altering 
its  passive  dynamics.  Through  consideration  of  the  coupling  between  the  structure  and  the 
adjacent  media,  this  provides  a  means  of  controlling  acoustic  sound  pressure  levels  or  adjacent 
flow  dynamics. 

A  great  deal  of  research  in  the  last  several  years  has  been  directed  toward  c^uestions  re¬ 
garding  the  modeling  of  piezoceramic  patch  interactions  with  underlying  structures  (see  [11] 
and  the  references  therein),  strategies  for  determination  of  optimal  patch  location  and  number, 
and  control  techniques  which  utilize  the  patches  as  sensors  and  actuators.  Due  to  the  steady 
state  periodic  nature  of  the  dynamics  in  many  structural  and  structural  acoustic  systems,  a 
large  number  of  the  current  control  methods  are  based  upon  freciuency  response  input/output 
analysis.  For  example,  Fuller  et  al.  have  employed  a  feedforward  filtered  X  version  of  an 
adaptive  LMS  algorithm  to  control  a  simply  supiDorted  plate  using  piezoceramic  actuators  [18]. 
While  such  techniciues  have  proven  cpiite  successful  for  controlling  steady  state  vibrations,  they 
do  not  have  the  capability  for  direct  control  of  transient  responses.  Other  successfully  im¬ 
plemented  methods  employing  piezoelectric  actuators  to  actively  control  structural  vibrations 
include  feedthrough  techniques  [17]  and  velocity  feedback  techniques  [1,  24].  In  general,  these 
methods  are  based  upon  modal  technicpies  and  are  designed  to  control  purely  steady  state 
responses.  An  exception  to  this  are  the  experimental  results  reported  in  [24]  in  which  tran¬ 
sient  plate  vibrations,  generated  by  an  impact  hammer,  were  reduced  using  an  analog  velocity 
feedback  circuit. 

Similar  studies  have  demonstrated  the  experimental  success  of  using  surface-mounted  piezo¬ 
ceramic  patches  to  reduce  structure-borne  noise  in  structural  acoustic  systems  [19,  20].  The 
emphasis  in  these  studies  was  again  on  using  frequency  input/output  analysis  to  control  steady 
state  dynamics. 
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An  alternative  approach  to  controlling  structural  vibrations  and  sound  pressure  levels  in 
structural  acoustic  systems  is  through  the  use  of  PDE-based  feedback  control  methods.  Analysis 
and  numerical  studies  demonstrating  these  technicpes  for  structural  applications  can  be  found 
in  [6,  7,  8,  12]  with  corresponding  results  for  structural  acoustic  systems  given  in  [4,  5,  10].  These 
techniques  start  with  an  infinite  dimensional  PDE  model  for  the  system  under  consideration. 
When  developing  such  models,  care  should  be  taken  to  incorporate  not  only  the  contributions 
due  to  the  piezoceramic  patches  but  also  dynamics  due  to  inexact  boundary  conditions  [13], 
coupling  with  adjacent  acoustic  fields  [5],  as  well  as  any  other  physical  phenomena  which  affect 
the  dynamics  of  the  structure.  In  this  setting,  mathematical  issues  such  as  model  well-posedness 
and  approximation  issues  concerning  simulations,  parameter  estimation  and  control  can  be 
addressed. 

By  approaching  the  problems  in  this  manner,  one  can  avoid  the  difficulties  caused  in  purely 
modal  methods  by  patch  contributions,  coupling  between  components,  and  inexact  boundary 
conditions.  Moreover,  by  combining  the  PDE  model  with  appropriate  time-dependent  feedback 
control  theory,  one  obtains  a  method  which  is  equally  applicable  for  controlling  transient  or 
steady  state  vibrations. 

In  this  paper,  the  experimental  implementation  of  such  a  PDE-based  control  method  is  con¬ 
sidered.  While  the  motivating  application  involves  the  control  of  vibration  levels  for  a  circular 
plate  through  the  excitation  of  surface-mounted  piezoceramic  patches,  the  general  techniques 
described  here  will  extend  to  a  variety  of  applications.  Following  a  brief  discussion  regarding 
the  model  and  a  Fourier/Galerkin  scheme  used  to  discretize  it,  relevant  feedback  control  theory 
is  discussed.  In  the  discussion  of  the  continuous  and  discrete  time  control  results,  two  cases  are 
considered;  namely,  the  control  of  plate  vibrations  in  the  absence  of  a  primary  input  force  and 
the  control  of  a  plate  driven  by  a  periodic  exogenous  force.  Implementation  issues  such  as  the 
effects  of  phase  shifts  and  delays  due  to  hardware  are  discussed  and  the  experimental  setup  is 
briefly  described.  Finally,  experimental  results  demonstrating  the  transient  and  steady  state 
control  results  are  presented.  These  demonstrate  the  effectiveness  of  the  PDE-based  controller 
for  this  system  and  indicate  the  potential  of  these  control  technicpies  for  reducing  transient  and 
steady  state  dynamics  in  other  structural  and  structural  acoustic  systems. 

Finally,  we  offer  comments  on  the  nature  of  our  contributions  here  to  the  literature.  This 
paper  does  not  contain  any  new  theoretical  results;  it  reports  on  our  successful  use  of  PDE- 
based  methods  in  experiments  at  NASA  Langley  Research  Center.  The  methods  (which  are 
theoretically  sound)  for  parameter  estimation  and  feedback  control  are  based  on  approximation 
theory  developed  (by  us  and  many  others  in  the  PDE  control  community)  during  the  past 
several  decades.  All  necessary  theorems  on  convergence  of  finite  dimensional  parameter  esti¬ 
mates,  gains,  filters,  controls,  observers,  tracking  variables,  et  cetera,  needed  for  the  examples 
treated  here  either  have  ai^peared  or  will  soon  appear  in  the  research  literature.  Many  of  these 
theoretical  results  on  PDE-based  control  have  been  largely  viewed  (especially  in  the  engineer¬ 
ing  community)  as  non-implementable  and,  hence,  as  somewhat  irrelevant  to  applied  scientists 
and  engineers.  The  present  manuscript  refutes  this  notion  and  validates  the  practicality  and 
importance  of  much  of  the  theoretical  efforts  over  the  past  years  on  PDE-based  control.  Our 
group  is  the  first,  to  our  knowledge,  to  provide  irrefutable  evidence  that  one  can  start  from 
basic  physical  laws,  derive  careful  infinite  dimensional  distributed  parameter  or  PDE  control 
models,  and  successfully  implement  PDE  control  methods  (with  the  approximations  necessary 
to  obtain  finite  dimensional  controls)  which  behave  as  theory  and  simulations  predict. 
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2  Circular  Plate  Model 


The  structure  under  consideration  is  a  thin  circular  plate  with  ^  sets  of  piezoceramic  patches 
bonded  to  the  plate  either  singly  or  in  pairs  as  depicted  in  Figure  2.1.  Throughout  this  discus¬ 
sion,  the  radius  and  thickness  of  the  plate  are  denoted  by  a  and  respectively.  Moreover,  the 
density.  Young’s  modulus,  Poisson  ratio,  Kelvin- Voigt  and  air  damping  parameters  are  given 
by  p,  E,  u,  Co  and  7  while  iv  and  g  are  used  to  denote  the  transverse  displacement  and  external 
force  on  the  plate.  Finally,  the  region  occupied  by  the  unstrained  neutral  surface  of  the  plate 
is  indicated  by  Fq. 

Equations  of  motion  for  the  plate  can  be  determined  from  both  Newtonian  (force  and 
moment  balancing)  and  Hamiltonian  (energy  formulation)  principles  and  we  summarize  both 
approaches  here.  The  presentation  will  be  for  a  general  Kirchhoff  plate  with  potentially  non- 
axiisymmetric  responses. 


Figure  2.1.  Thin  circular  plate  with  surface-mounted  piezoceramic  patches. 


2.1  Strong  Form  of  Plate  Model 

Considering  first  the  model  which  derives  from  Newtonian  principles,  we  let  M,.,  Mg,  M^e  denote 
internal  moments  and  {Mr)pe,  denote  external  moments  generated  by  the  piezoceramic 

patches.  As  detailed  in  [9,  29],  for  a  structure  with  patch  pairs  activated  out-of-phase,  the 


internal  moments  are  given  by 


Mr  =  -D 


Mf,  =  -D 


d^W  V  dw  V  d^W 

dr‘2  r  dr  86^ 


1  dio  1  d^i 


( 1  d^w 

Mre  =  -D{l-v)[-^, 


1  dw 


/  d^w  y  d^w  V  d^w  \ 
'^\drm^^dm'^^'d¥m) 

( 1  d^w  1  d^w  d^vj  '' 
r  ^  ^  dQm  ^  drm, 

(i  d\v  1  d^io 


where  D  =  i2(T-^)  denotes  the  plate  stiffness.  The  external  moments  generated  by  the  s  patch 
pairs  have  the  form 


iMr)pe  =  iMe)pe  =  -  ^  ICf  Hi{t)xi{r,  0)  ,  Xi(r,  6) 


1  ,  (r,  patch  region 

0  ,  elsewhere 


where  Ui{t)  is  the  voltage  into  the  patch  and  /Cf  is  a  parameter  which  depends  on  the  geom¬ 
etry,  piezoceramic  material  properties  and  piezoelectric  strain  constants  (see  [11]  for  details). 
The  internal  and  external  moments  can  then  be  combined  to  yield  the  general  plate  moments 

Mr  =  Mr-{Mr)pe' 

Me  =  Me  —  {Me)pe 
Mre  —  Mre  • 

For  a.  clamped  plate  with  initial  displacement  woir,$)  and  velocitj^  ivi{r,0),  force  and  mo¬ 
ment  balancing  yield  the  equations 


d^w  dw 

'''■ac 


d^Mr  2dMr  ^  IdMe 

dp2  p  Qp  p  Qp 


'  dt  dr'^  r  dr  r 
2  d‘^Mre  _  2  dMre  _  1  d'^Me 
r  drd0  7’^  d0  d^^ 


=  gii-,i'-,0) 


0  <  6'  <  2  7r 
0  <  r  <  a 


ui(t,  a,  0)  =  a,0)  =  0 

dr 


|r(j(0,r,^)  =  u;o(r,^>)  ,  ^(0,  r,  6»)  =  u)i(r,  6>)  . 

as  the  strong  form  of  the  plate  model  (see  [9,  27]  for  details).  From  an  apjslications  perspective, 
the  following  observations  can  be  made  regarding  this  model. 

1.  It  is  first  noted  that  in  many  applications,  it  is  nearly  impossible  to  maintain  the  truly  fixed 
(zero  displacement  and  slope)  boundary  conditions  specified  in  the  model  (2.1).  For  the 
experimental  plate  which  we  used,  parameter  estimation  results  indicated  minimal  energy 
loss  through  the  boundary  conditions  and  an  adec^uate  fit  of  the  model  (2.1)  was  obtained 
(see  the  identification  results  in  [2,  3]).  For  plates  in  which  the  boundary  clamping  is  less 
secure,  a.  model  for  imperfectly  clamped  boundaries  such  as  that  presented  in  [13]  should 
be  used  in  order  to  obtain  a  model  fit  which  is  adeciuate  for  control  applications. 
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2.  As  discussed  in  [11],  the  plate  parameters  p,D,u  and  cd  are  discontinuous  due  to  the 
presence  and  differing  material  properties  of  the  patches.  Moreover,  while  ‘handbook’ 
values  can  be  determined  for  p,  D  and  n  for  a  plate  which  is  devoid  of  patches,  those  values 
usually  cannot  be  used  in  the  final  system  model  with  any  accuracy  due  to  nonuniformities 
in  the  plate  and  boundary  conditions,  variations  in  materials,  and  contributions  due  to 
the  pre.sence  of  the  i^atches.  In  applications,  such  as  that  in  this  paper,  these  parameters 
are  estimated  using  fit-to-data  techniques  (see  [2,  3]  for  the  parameter  estimation  results 
pertaining  to  the  plate  used  here). 

3.  The  input  parameters  /Cf  are  discontinuous  since  they  are  nonzero  only  over  the  regions 
of  the  patches.  While  expressions  for  these  constants  are  derived  in  [11],  they  too  must 
be  estimated  in  applications  due  to  manufacturing  variations  in  the  patches. 

It  is  readily  noted  in  the  strong  form  of  the  modeling  equations  that  the  discontinuous  plate 
parameters  and  patch  input  terms  are  differentiated,  thus  leading  to  derivatives  of  the  Dirac  6 
‘function’.  The  difficulties  associated  with  this  formulation  are  avoided  in  the  weak  form  of  the 
modeling  ec|uations  which  is  presented  in  the  next  section. 


2.2  Weak  Form  of  Plate  Model 

To  provide  a  framework  which  facilitates  analysis,  approximation  and  implementation,  it  is 
advantageous  to  consider  a  weak  form  of  the  modeling  ecpiations.  Such  a  formulation  can  be 
determined  directly  from  Hamiltonian  (energy)  principles  and  is  equivalent  to  that  obtained  by 
integration  by  parts  after  multiplication  of  the  strong  form  by  suitably  smooth  test  functions. 

For  the  plate  problem  under  consideration,  a  suitable  space  V  of  test  functions  is  the  subset 
of  the  Sobolev  space  H^{To)  which  satisfies  the  essential  boundary  conditions  =  ^  =  0  at 
r  =  a.  As  detailed  in  [9,  29],  a  weak  or  variational  form  of  the  equations  of  motion  for  the  plate 
is 


f  ,  f  f  ,,  ,  f  1  -  ,  r  ,  d'^1]  , 


2/  -At 

Jto 


e 


d’^T]  dr] 

drdO  do 


^/Cfiq(t)xi(r,6>)V2j;d7  +  /  gi]d') 

•  1 


for  all  rj  G  V .  The  overbar  here  denotes  complex  conjugation  and  the  differential  is  dy  =  rdOdr. 
It  is  easily  noted  that  in  this  form,  derivatives  that  were  originally  applied  to  moments  have  been 
transferred  onto  test  functions.  This  eliminates  the  difficulties  associated  with  differentiating 
the  piecewise  constant  parameters  p,  D,  v  and  cd  found  in  the  internal  moments  as  well  as  the 
discontinuous  injaut  parameters  /Cf ,  •  •  • ,  /Cf . 


2.3  State  Approximation 

As  discussed  in  [9,  29],  an  appropriate  choice  for  the  basis  and  Fourier-Galerkin  expansion 
of  the  plate  displacement,  when  considering  clamped  boundary  conditions,  is  B^(r,0)  = 


M  AT 

i:  E’‘’l('H"‘‘'K(>'K"'’  =  E<‘>f(i)B^(''.»)-  (2-3) 

m=—M  71—1  k—1 

Here  5”^?;)  is  the  modified  cubic  spline  satisfying  B^(a)  =  =  0  with  the  condition 

=  0  being  enforced  when  m.  =  0  (this  latter  condition  guarantees  differentiability  at  the 
origin  and  implies  that 

f  N  ,  m  =  0 
=  < 

\  +  1  ,  ?77.  7^  0 

where  N  denotes  the  number  of  modified  cubic  splines).  The  total  number  of  plate  basis 
functions  is  Af  =  {2M  +  1)(A^  +  1)  -  1.  As  discussed  in  the  [9,  29],  the  inclusion  of  the  weighting 
term  with 

[  0  ,  ?77.  =  0 

"  I  1  ,  777  7^  0 

is  motivated  by  the  asymptotic  behavior  of  the  Bessel  functions,  which  make  up  the  analytic 
plate  solution,  as  r  — >•  0.  It  also  serves  to  ensure  the  uniqueness  of  the  solution  at  the  origin. 
The  Fourier  coefficient  in  the  weight  is  truncated  to  control  the  conditioning  of  the  mass  and 
stiffness  matrices  (see  the  examples  in  [9]). 

To  obtain  a  matrix  system  (again,  see  [9,  29]  for  a  careful  derivation  with  complete  details), 
the  Af  dimensional  approximating  subspace  is  taken  to  be  =  span{B'^}  and  the  product 
space  for  the  usual  corresi^onding  first-order  vector  system  is  x  .  The  restriction  of  the 
first-order  form  for  the  infinite- dimensional  system  (2.2)  to  the  space  x  then  yields  the 
matrix  equation 


■  0  j  r  1  r  0  j  0  r  0 

0  .  -K^  .  _  ^^{t)  _  _  .  1. 

'  0  ]  r  77^(0)  1  \  ' 

.  0  \  7?^(0)  ~  [g^  . 


where  d^{t)  =  [u7f^(t ),  (t),  •  •  • ,  denotes  the  column  Af  vector  containing  I  he  approx¬ 

imate  state  coefficients  (see  (2.3)).  The  component  matrices  and  vectors  are  given  by 


—  Kdi  +  +  a  +  Ada  +  Ad^ 

=  KcnA  +  +  A',^3  +  A'e^4  + 

JFo 

=7  phBif^di, 

>•  J  AA;  7ro 

Kl, =(”’«’<)>•  ’  KL =(«'»’<>« 
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where 


[Knikk  =  I 

JTo 

(A-Ddy-  = 


D 


D 


D 


d^B{f  V  dBt'  V  d'^Bf 


dr 


r>2 


+ 


r  dr 


+ 


do^ 


d^Bf 

dr'^ 


dl  , 


1  dB^  ,  1  d‘^B^  .  vd'^Bif 


+ 


+ 


J.2  Q^>  j,3  QQ2  Qj.2 

1  dB^  1  d'^B^  V  d^Bf 


dBf 

<9r 


d-f  , 


d^Bf  ^ 


[A'mlf,/;  =  2  /  D{l-v) 
Jtq 

[^'l?4fc  =  2/  DH-v) 

JTq 


1  d'^-B^  1  dBf 


drdd  do 


1  d^Bf  1 


7*3 


dr  do 


+ 


"4  c>0 


drdO 
dBf 

do 


dl  , 


d'y  . 


The  index  ranges  here  are  k^  £  =  1,  *  •  * ,  A^-  The  matrices  Ac^i  —  defined  similarly  with 

the  inclusion  of  the  parameter  cjj  in  the  various  integrals.  Finally,  we  remind  the  reader  that 
p,  D,  If  and  are  piecewise  constant  in  these  definitions  due  to  the  presence  of  the  patches. 

For  application  purpose,  it  useful  to  note  that  the  matrix  system  for  the  plate  can  thus  be 
written  as  the  Cauchy  system 


y^{t)  =  +  B"u{t)  +  y"(t) 

y^iO)  =  yo 

where  y^{t)  =  ■  ■  ,wj![r{t),ibf{t),- •  ■  ,wfrit)f  denotes  the  column 

2^/"  vector  containing  the  generalized  Fourier  coefficients  for  the  approximate  displacement  and 
velocity.  In  this  form,  the  control  problem  can  be  i-eadily  discussed. 


3  Continuous  Time  Control  Problem 

In  the  preceding  discussion  leading  from  the  infinite  dimensional  PDF  model  to  the  finite  dimen¬ 
sional  matrix  approximation,  the  superscript  V  was  used  to  denote  the  level  of  discretization; 
i.e.,  the  number  of  Fourier/spline  basis  elements  used  to  approximate  the  state.  This  notation 
is  standard  in  the  theory  of  finite  element  and  spline  approximations  of  infinite  dimensional 
systems.  In  finite  dimensional  control  theory,  however,  the  level  of  discretization  is  typically 
fixed  and  these  subscript  V’s  are  usually  omitted  to  simplify  notation.  We  will  do  the  same  in 
this  and  subsecpent  sections  so  as  to  remain  consistent  with  standard  control  notation. 

3.1  Initial  Displacement  and  Velocity  —  LQG  Control  Law 

We  consider  first  the  V-dimensional  systems 

y{t)  =  Ay{t)  -h  Bu(t}  ,  y{0)  =  yo 

yob(t)=Cy{t)  (3.1) 

zit)  =  Hy[t)  +  Gti{t) 
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where  yob  denotes  observations  in  ]R^  and  C  \s  a.  P  X  Af  observation  matrix  whose  structure 
is  determined  by  the  manner  and  number  of  observations  being  used  (the  specific  C  matrix 
used  in  the  plate  experiments  is  described  in  Section  4.2).  Moreover,  ^  6  IR^  denotes  the 
performance  output  obtained  under  the  assumption  that  G  and  H  are  time  invariant  matrices 
satisfying  H^G  =  0.  In  the  event  that  P  =  Af  and  G  is  an  identity,  the  optimal  control  u  can  be 
obtained  from  standard  linear  cjuadratic  regulator  (LQR)  optimal  control  theory.  The  number 
of  observations  P  is  usually  limited,  however,  and  we  concentrate  instead  on  the  case  P  <  Af 
which  occurs  when  the  full  state  is  unavailable  and  must  be  reconstructed  using  a  compensator 
(e.g.,  see  [25]). 

The  general  control  problem  for  this  case  consists  of  determining  the  voltage  u  which  mini¬ 
mizes  the  performance  index  (or  cost  functional) 

J[u)  =  f  \zit)f  dt 

Jo 

/•OO 

=  /  {{Qy(t),yit))  + {Ru{t),u{t))}dt 

Jo 

subject  to  (3.1).  The  A'*  x  A/”  matrix  Q  can  be  chosen  to  satisfy  various  design  criteria  including 
frequency  windowing,  the  weighting  of  various  state  components,  or  minimization  of  certain 
energy  measures.  The  Q  matrix  used  here  was  chosen  using  energy  considerations  and  con¬ 
struction  details  can  be  found  in  Section  4.2.  The  s  x  s  matrix  R  weights  the  voltages  to  the 
various  patches  or  patch  pairs. 

Because  full  state  information  is  not  available  in  most  applications,  the  state  must  be  esti¬ 
mated  or  reconstructed  from  observations  before  a  controlling  voltage  can  be  determined.  We 
consider  here  a  full  order  compensator  or  observer  of  Luenberger  type  [25]  and  refer  the  reader 
to  [16,  22]  for  details  on  reduced  order  observers. 

The  comiDensator  or  reconstructed  state  satisfies  the  matrix  system 

ydO  =  Ayc(t)  Bu{t)  +  F  [yob{t)  —  Cydt)] 

ydO)  =  ?/co 


with  the  optimal  voltage 


u{t)  =  -Kydt) 


where  F  and  K  denote  the  compensator  and  feedback  gains,  respectively.  We 
and  F  are  chosen  so  that  the  reconstruction  error  \y{t)  —  ?/c(A)l  0  as  t  ^  oo. 

observability  and  controllability  hypotheses  (see  [25]),  the  optimal  feedback  and 
gains  are  given  bv 

K  = 

F  =  PC^R-^ 


note  that  K 
Under  usual 
compensator 


(3.2) 


where  11  and  P  are  unique  nonnegative-definite  solutions  to  the  following  feedback  (regulator) 
and  compensator  (observer)  algebraic  Riccati  eciuations 


IIA  +  A^U  -  UBR-^B^n  +  Q  =  0 
PA^  A  AP-  PG'^R-^CP  +  (?  =  0  , 
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respectively.  As  was  the  case  with  the  matrices  Q  and  i?,  the  matrices  Q  and  R  are  design 
criteria  for  the  specific  control  application  under  consideration  (specific  choices  used  in  the  plate 
experiments  are  summarized  in  Section  4.2).  We  point  out  that  in  terms  of  the  component 
matrices  and  control  voltage,  the  compensator  can  be  expressed  as 


ydt) 


A  -  BR-^B^n]  yc(t)  +  PC^R-^C  [y(t)  -  yM  . 


The  control  law  just  described  must  be  implemented  in  real  time  in  order  to  be  a  viable 
method  for  reducing  vibrations  in  physical  structures.  To  facilitate  implementation,  it  is  pru¬ 
dent  to  calculate  offline  as  many  components  as  possible  and  then  treat  those  pre-calculated 
components  as  filters  when  performing  online  computations.  The  method  for  continuous  time 
implementation  is  summarized  and  offline  and  online  components  are  categorized  in  the  follow¬ 
ing  algorithm. 


Algorithm  3.1.1:  Continuous  Time  Control  of  Initial  Displacement 

and  Velocity 


Offline  (f) 

(fi) 

(Hi) 

Construct  matrices  A,  fi,  C,  Q,  R,  R 

Solve  Riccati  equations  (3.3)  for  II  and  P 
Construct  K  =  R~^B^Il  and  F  =  PC^R~^ 

Online  (?) 

Collect  data  yob{t)  =  Cy{t) 

(ii) 

Solve  the  ODE  system 

ydt)  =  [A-BK  -FC]ydt)  +  Fyodt) 

=  Acydt)  +  Fyobit) 

{Hi) 

Calculate  the  voltage  n{t)  =  —Kydt) 

We  note  that  the  expensive  (time  consuming)  calculation  of  the  component  matrices  A,  B,Q^R^ 
C,  Q,  R  and  Riccati  solutions  If,  P  is  performed  offline  with  the  results  loaded  into  the  control 
code  as  datafiles.  This  leaves  the  integration  of  the  system  ydt)  =  Acydi)  +  Pyob{t)  as  the 
primary  computation  to  be  performed  during  implementation.  Issues  regarding  the  numerical 
integration  of  the  system  as  well  as  the  effects  of  discrete  time  calculations  will  be  discussed  in 
Section  4,  and  a  discrete  version  of  Algorithm  3.1.1  is  summarized  in  Algorithm  4.1.1.  Finally, 
the  observed  system,  performance  index  and  state  estimator  and  LQG  control  law  for  the  system 
are  summarized  in  Algorithm  3.4.1  in  Section  3.4. 

While  this  compensator  does  i^rovide  the  desired  performance,  it  may  lack  robustness  in 
some  applications.  In  cases  where  added  robustness  with  regard  to  certain  types  of  system  or 
observation  noise  and  modeling  errors  is  recpiired,  an  if^/MinMax  compensator  of  the  type 
described  in  the  next  section  can  be  used. 
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3.2  Initial  Displacement  and  Velocity  —  iJ°^/MinMax  Control  Law 

We  consider  now  the  design  of  a  dynamic  compensator  which  is  robust  with  respect  to  certain 
types  of  state  and  measurement  uncertainties  or  disturbances  (see  [15]  for  details).  To  incorpo¬ 
rate  such  uncertainties,  we  let  w{t)  G  IR^  denote  input  and  output  disturbances.  The  system, 
with  no  exogenous  force,  is  then  given  by 

y{t)  =  Ay{t)  +  Bu{t)  -f  Du){t) 

yob{t)  =  Cy{t)  -h  Ew{t) 

z[t)  =  Hy{t)  +  Gu{t)  . 

For  this  discussion,  we  assume  that  the  input  and  output  disturbance  are  independent  and 
hence  DE^  =  0.  This  is  a  matter  of  convenience  and  the  dependent  case  can  be  handled 
similarly  after  slight  modihcations  are  made  (see  [15]). 

In  this  case,  the  min/max  optimization  problem  leading  to  the  controller  consists  of  finding 
a  controller  u*  G  U  =  T^(0,  oo;  IR,*)  and  disturbance  w*  G  IT  =  T^(0,  oo;lR®)  such  that 

J*  —  inf  sup  J^{u,w)  =  J^{u*,w*) 

^  u^U 

for  the  disturbance  augmented  functional 

J^{u,w)  =  ^  \w{t)f]  dt 

=  {{Qy(t),y{t))  + ^ 

As  noted  in  [12,  15],  the  results  from  this  optimization  problem  yield  a  bound  7  for  the 
norm  of  the  transfer  function  from  disturbance  C{w)  to  the  performance  output  C{z)  wdrere  C 
denotes  the  Laplace  transform. 

Under  the  assumption  that  the  i^air  {A,B)  is  stabilizable,  {A,C)  is  detectable,  (A,  G)  is 
controllable,  and  (A,  H)  is  observable,  one  can  prove  the  existence  of  (minimal)  positive  definite 
solutions  n  and  P  to  the  algebraic  Riccati  equations 

RA  -f  A^n  -  n  BR-^B^  -  7-^Q]  n  +  G  =  0 

PA^  +  AP-P  [c’^ir^C  -  7-2g]  p  +  g  =  0 

for  a  given  attenuation  7  >  0.  Moreover,  if  the  spectral  radius  p  of  PII  satisfies  the  condition 

p(Pn)  <  7^  or  n  -  q'^p-'  <  0, 

then  there  exists  a  unique  optimal  controller 

u*{t)  =  -R-^B'^B.ydt)  . 

The  state  estimator  ydi)  £  IR"^  satisfies 

ydi)  =  Acycit)  +  Fyoh{t) 

?/c(0)  =  yco 
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where 


A,  =  A- BK-  FC  +  7“^ On 


F  = 


I  -  j-^PU  "  PC^  R 


The  implementation  issues  concerning  the  method  are  similar  to  those  discussed  in  the 
algorithm  for  the  LQG  controller  but  determination  of  the  Riccati  solutions  If,  P  and  filter  F 
are  complicated  by  the  fact  that  a  suitable  design  value  of  7  must  be  determined  before  matrix 
calculations  can  proceed.  Fortunately,  these  calculations  can  be  performed  offline  and  resulting 
matrices  input  as  data  files  for  the  online  computations.  Hence  the  actual  online  controller  can 
run  at  the  same  rate  as  that  obtained  using  the  LQG  methodology. 

The  form  of  the  observed  system,  the  extended  performance  index  and  the  F/'^/MinMax 
control  law  are  summarized  in  Algorithm  3.4.2  in  Section  3.4.  The  various  components  can  be 
compared  with  the  LQG  components  summarized  previously  as  well  as  those  which  arise  when 
a  periodic  exogenous  force  is  considered. 


3.3  Periodic  Primary  Excitation 

For  the  case  in  which  a  periodic  exogenous  force  drives  the  system,  knowledge  of  that  force 
can  be  used  extend  previously  discussed  results  to  include  the  effects  of  periodicity.  The  Af- 
dimensional  system  in  this  case  is 

y(t)  =  Ayit)-\- git)  ,  j/(0)  =  y(T) 

(3.4) 

yob(t)  =  Cyit) 


where  g(t)  G  IR’^  is  periodic  with  period  r.  That  periodicity  is  then  reflected  in  the  performance 
index 

=  /  {{Qyit),y{t))  +  {R'it{i)^u(t))}  dt 

Jo 

which  is  minimized  subject  to  (3.4). 

With  K  and  F  defined  in  (3.2),  the  reconstructed  state  in  this  case  satisfies  the  system 
ydt)  =  Aydt)  +  Bu{t)  +  F  [yobit)  -  Cydt)]  +  g(t) 

(3.5) 

ydO)  =  ?/c(r) 

with  the  optimal  voltage  given  by 

u(t)  =  -Kydt)  +  R~^B'^r{t)  .  (3.6) 


Here  r  is  a  tracking  variable  defined  by  the  system 

r(i)  =  -[.4-BA']’'r(()  +  n9(() 
r(0)  =  r(r) 

where  H  solves  the  first  of  the  algebraic  Riccati  ecpiations  (3.3).  We  point  out  that  in  this 
case,  the  voltage  contains  two  contributions.  The  first  incorporates  transient  information  bj^ 
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feeding  back  state  estimates  while  the  tracking  component  incorporates  information  regarding 
the  periodic  force. 

Combining  (3.5)  and  (3.6)  yields  the  single  expression 


il,{t)  =  [A-  BK]  ydt)  +  FC  [y{t)  -  y^t)]  +  BR  ^B'^r{t)  +  g{t) 

2/c(0)  =  ydr)  . 


(3.8) 


for  the  state  estimator  in  terms  of  the  tracking  variable.  The  state  estimate  at  time  t  is  then 
obtained  by  integrating  (3.7)  and  (3.8)  after  solving  the  necessary  Riccati  equations. 

As  was  the  case  when  considering  control  of  the  unforced  system  with  initial  displacement 
and  velocity,  the  computations  can  be  categorized  with  respect  to  those  which  can  be  performed 
offline  and  those  which  must  be  done  online.  The  following  algorithm  summarizes  the  continuous 
time  control  method  for  periodic  excitation  and  categorizes  the  offline  and  online  components. 


Algorithm  3.3.1:  Continuous  Time  Control  with  Periodic 

Exogenous  Force 


Offline  (i) 

(li) 

{Hi) 

Construct  matrices  A,  B,  C,  Q,  R,  Q,  R 

Solve  Riccati  equations  (3.3)  for  11  and  P 

Construct  K  —  R~^B^Il  and  F  =  PC^R"^ 

Online  (i) 

Collect  data  yob{t)  —  Cy{t)  and  force  measurements  g{t) 

(?:{) 

Solve  the  ODE  systems 

r(t)  =  -[A  -  BKf  r{t)  +  Ug{t) 

=  Atrrit)+g{t) 

ydt)  =  [A~BK-FC]ydt)  +  Fyodt)  +  BR-^B'^r{t)  +  git) 

=  Acydt)  +  Fyob{t)  A  Arr{t)  +  g{t) 

(iti) 

Calculate  the  voltage  u{t)  =  —Kydt)  +  R~^B^r{t) 

In  this  case,  the  tracking  equation  r{t)  =  Atrr{t)  +  g{t)  must  be  solved  before  the  state  can 
be  estimated  by  integrating  the  system  ydt)  =  Acydt)  +  Fyob{t)  +  A,.r(f)  +  g(t)-  Details 
regarding  the  solution  of  these  two  systems  with  discrete  data  and  force  measurements  are 
given  in  Algorithm  4.1.2  in  the  next  section. 

Finally,  the  governing  equations  and  component  matrices  for  the  LQG  controller  are  sum¬ 
marized  in  Algorithm  3.4.3  of  Section  3.4.  Expressions  similar  to  those  in  Section  3.2  arise 
when  an  compensator  is  considered  for  the  problem  and  the  resulting  control  law 

is  summarized  in  Algorithm  3.4.4.  We  also  refer  the  reader  to  [4]  where  details  concerning  the 
design  of  an  If^'/MiiiMax  compensator  for  a  structural  acoustic  system  that  is  subjected  to  a 
periodic  exogenous  force  is  considered. 
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3.4  Summary  of  Systems,  Performance  Indices  and  Control  Laws 

The  modeling  state  equations,  observed  system,  performance  index  and  compensator  design 
for  the  previouslj^  discussed  cases  are  summarized  in  this  section.  Specifically,  the  LQG  and 
i/'^/MinMax  control  laws  are  summarized  for  the  system  with  no  primary  force  as  well  as 
the  system  subjected  to  a  periodic  exogenous  force.  Details  regarding  the  construction  of  the 
individual  components  can  be  found  in  the  preceding  section  or  in  the  implementation  section 
which  follows. 


Algorithm  3.4.1:  LQG  Control  Law  for  Unforced  System 


(1)  Observed  System 

y{t)  =  Ay{t)  +  Bu{t) 
yoh{i)  =  Cy{t) 
z{t)  =  Hy(t)  +  Gu(t) 

y  e  «  6  E* 

yob  e  E'’ 

zeW  ,  h'^G  =  0 

(2)  Performance  Index 

./(«)=  [  \z(t)fdt 

Jo 

=  [  {{Qy(t),y{t))  + {Ru{t),u{t))}dt 

Jo 

Q  =  H'^H  ,  R  =  G^G 

(3)  State  Estimator  and  Control  Law 

ydO  =  Acycit)  +  Fyohit) 

Ac  =  A- FC-  BK 

K  = 

F  =  PC'^R-^ 

iLA  +  -  uBR-^B^n  +  g  =  0 

PA^  +  AP-  PC^R-^CP  +  Q  =  0 

Controlling  Voltage 

u(t)  =  -Kydt) 
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Algorithm  3.4.2:  iJ°^/MinMax  Control  Law  for  Unforced  System 


Algorithm  3.4.3:  LQG  Control  Law  for  Periodic  Exogenous  Force 

(1)  Observed  System 

y{t)  =  Ay{t)  +  Bu{t)  +  g{t) 

y(0)  =  y{T) 

yoh{t)  =  Cy{t) 

z{t)  =  Hy{t)^Gu{t) 

(2)  Performance  Index 

J(u)=  [\z{t)\Ut 

JO 

=  /  {{QyiOtViO)  ^  Q  =  H^H  ,  R=G^G 

Jo 

(3)  State  Estimator  and  Control  Law 

ydt)  =  AcVcit)  +  Fyobit)  +  A,r(t)  +  g{t) 

Ac  =  A- FC  -  BK 
Ar  =  BR-^B^ 

K  =  R-^B'^B 
F  =  PC'^R-^ 

n.4  +  A^B  -  BBR-^B'^B  +  Q  =  0 

PA^  +  AP  -  PC^R-^CP  +  Q  =  0 

Tracking  Equation 

r[t)  =  Atrr(t)  P  cj{t) 

Atr  =  -[A-  Bhf 
g{t)  =  n^(t) 

Control  Voltage 

u{t)  —  —Kydt)  +  R~^  B^r{t) 


y  e  u  e  E" 
yob  e  E'^ 

^  6  E^  H'^G  =  0 


Algorithm  3.4.4:  iJ°^/MinMax  Control  Law  for  Periodic  Exogenous 

Force 

(1)  Observed  System 

y{t)  =  Ay{t)  +  Bu{t)  +  Dw{t)  +  g{t)  y  6  IR^^,  u  G  ]R®,  tw  G  IR® 

3/(0)  =  3/(r) 

yob{t)  =  Cy{t)  +  Ew{t)  yob  G  IR.^  ,  DE'^  =  0 

z{t)  =  Hy(t)  +  Gu{t)  2  G  E’’  ,  H'^G  =  0 

Q  =  >  0  ,  R  =  G^G  >  0 

Q  =  DD^  >  0  ,  R  =  EE^  >  0 

(2)  Performance  Index 

j^(u,ic)  =  dt 

(3)  State  Estimator  and  Control  Law 

ycit)  =  Acydt)  +  Eyob{t)  +  Arr{t)  +  g{t) 

Ao  =  A-  EC  -  BK  + 

Ar  =  BR-^B^  -  ^-'^Q 

K  =  R~^B^U 

F  =[i PC^R-^ 

EA  +  A'^U  -  n  BR-^B^  -  7-"g]  R  +  Q  =  0 
PA^ +  AP-P[C^R-^C--f-^Q]p  +  Q  =  0 

Tracking  Equation 

r{t)  =  Atrr(t)  +  g(t) 

Atr  =  -  [A-BK  -^-^QeY 

g{t)  =  Eg{t) 

Control  Voltage 

u{t)  =  —Kydt)  +  R~^B^r(t) 


Note:  piPE)  <  or  H  -  <  0 


4  Discrete  Time  Control  Problem 


The  control  laws  summarized  in  the  last  section  were  derived  under  the  assumption  of  continuous 
time  sampling  of  observed  state  and  force  data.  Hence  it  was  assumed  that  yob{t)  and  g{t)  were 
available  for  all  t  within  the  temporal  interval  of  interest.  Moreover,  it  was  assumed  that 
y{t),yc(t)  and  r{t)  could  be  obtained  through  exact  integration  of  the  state,  state  estimator 
and  tracking  equations.  When  implementing  the  method,  however,  one  has  available  only 
discrete  data  values  and  the  differential  estimator  and  tracking  equations  must  be  numerically 
approximated  at  discrete  time  values.  The  manner  in  which  the  control  laws  are  implemented 
in  discrete  time  and  the  influence  of  this  discretization  on  the  overall  performance  are  discussed 
in  this  section. 

In  the  last  section,  it  was  demonstrated  that  in  the  case  with  no  primary  exogenous  force, 
the  state  equation,  state  estimator  and  controlling  voltage  for  the  LQG  and  if'^'/MinMax 
problems  had  the  form  summarized  in  column  1  of  Table  4.1.  The  corresponding  quantities  for 
a  system  subjected  to  a  periodic  exogenous  force  are  summarized  in  the  second  column.  Details 
regarding  the  component  matrices  can  be  found  in  the  previous  section. 


No  Primary  Input 

Periodic  Exogenous  Force 

State 

Equation 

y{t)  =  Ay(t)  +  Bu{.t)  ,  y(0)  =  yo 

Vobit)  =  Cy(t) 

y(t)  =  Ay{t)  +  Bu{t)  +  (j(t)  ,  t/fO)  =  2/(r) 

yobit)  =  Cy{t) 

State 

Estimator 

ydt)  =  Acydt)  +  Fyobit) 

2/c(0)  =  2/co 

ydt)  =  Acydt)  +  Fyob(t)  +  Arr{t)  +  g{t) 

2/c(0)  =  ydr) 

Tracking 

Equation 

f[t)  =  Atrr(t)  +  g(t) 

7'(0)  =  r(r) 

Control 

Voltage 

u{t)  =  -R~'^ 

u(i)  =  ~R~^ 

Table  4.1.  State,  state  estimator  and  controlling  voltage  for  the  systems  with  no  primary 
injrut  and  a  periodic  exogenous  force. 

The  relationship  of  these  components  to  the  experimental  plate  setup,  driven  by  a  periodic 
exogenous  force,  is  illustrated  in  Figure  4.1.  For  that  setup,  the  physical  process  consists  of  the 
clamped  circular  plate  with  attached  piezoceramic  patches.  The  state  equations  represent  the 
spatial  discretization  of  a  PDF  model  of  the  process  with  physical  parameters  estimated  using 
experimental  data  so  that  the  model  accurately  captures  the  plate  dynamics. 
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Physical  Process 


on 


Exogenous 

Force 


Observation  a(t) 


Exogenous  Force  g(t) 


Feedback 

u(t)  =  -Kyc(t)  +  R^B^r(t) 


State  Equations 
(Model  Process) 

y(t)  =  Ay(t)  +  Bu(t)  +  g(t) 

yob(t)  =  Cy(t) 


Filters 
DSP  Board 


PC  Computer 

(1)  Integrate  Acceleration  Data 
a(t)  ->  v(t) 

(2)  Determine  Tracking  Component 
;(t)  =  Atrr(t)+ng(t) 

(3)  Estimate  Process  from  Observations 

yc  (t)  =  Acyc  (t)  +  Fv(t)  +  br‘  B^ r(t)  +  g(t) 

(4)  Calculate  Voltage 

u(t)  =  -Kyc(t)  +  R’BVt) 


Figure  4.1.  Experimental  plate  setup,  modeling  equations  and  computations  necessary  for 
determination  of  controlling  voltage. 

We  point  out  that  when  modeling  the  process  or  performing  simulations,  the  observations 
yobii)  are  obtained  from  the  approximate  state  through  the  relation  yob{^)  =  Cy(t)  where  C 
is  the  matrix  observation  operator.  When  experimentally  implementing  the  method,  temporal 
data  from  the  process  is  used  to  determine  the  observations  yob{t)-  This  often  requires  processing 
either  through  digital  or  analog  filters  or  in  the  software.  In  the  circular  plate  experiments, 
accelerometer  data  a{t)  was  integrated  to  yield  velocity  measurements  v{t)  to  be  used  for 
the  observations  yob{t)  (see  the  final  subsection  of  this  section  for  details).  Once  yob[^‘)  has 
been  obtained,  the  state  is  estimated  by  numerically  integrating  the  estimator  equation  and  a 
controlling  voltage  is  calculated.  This  voltage  is  then  fed  back  to  the  process. 

Because  data  is  processed  in  a  digital  manner,  the  observations  yob  and  force  measurements 
g  are  obtained  only  at  discrete  times  tj.  The  rate  at  which  this  data,  can  be  sampled  is  governed 
by  the  data  acquisition  system,  the  software  being  used  and  the  number  of  calculations  required 
between  samples.  In  particular,  the  tracking  component  r(tj),  state  estimate  ydij)  and  voltage 
u{tj)  must  be  calculated  before  the  arrival  of  data  at  time  tj+i-  While  details  regarding  these 
calculations  are  postjDoned  until  Section  4.1,  it  is  noted  here  that  the  system  sizes  must  be 
minimized  in  order  to  permit  real-time  approximation  of  the  estimation  equation.  This  is 
a  major  motivation  for  choosing  ai^propriate,  accurate  approximation  techniques  to  siDatially 
discretize  the  modeling  PDE’s. 

As  indicated  previously,  the  open  and  closed  loop  process  dynamics  can  be  simulated  by 
using  the  state  equations  to  model  the  plate.  The  state  estimate  and  state  are  then  approxi- 
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mated  by  simultaneously  integrating  the  corresponding  equations.  The  integration  of  the  state 
equation  can  be  performed  using  any  sufficiently  accurate  ODE  routine  which  is  efficient  for  the 
system  under  consideration;  for  example,  a  variable  order,  variable  stepsize  method  was  used  to 
solve  the  stiff  system  which  arose  when  simulating  plate  dynamics.  If  simulations  demonstrating 
the  levels  of  control  that  can  be  obtained  under  “optimal”  conditions  are  desired,  the  estimator 
equation  can  be  integrated  nsing  the  same  high-order  routine.  On  the  other  hand,  simulations 
representing  the  attenuation  levels  that  can  be  expected  under  “implementation”  conditions  can 
be  obtained  by  incorporating  values  of  yob  and  g  calculated  at  discrete  times  and  approximating 
ydtj)  using  the  techniques  employed  when  implementing  the  method.  Simnlation  results  using 
both  techniques  can  be  found  in  [14]. 

4.1  Approximation  of  the  Estimator  and  Tracking  Equations 

In  order  to  obtain  tracking  values  and  state  estimates  to  be  used  when  calculating  controlling 
voltages,  the  solution  to  the  tracking  and  state  estimator  equations  must  be  numerically  ap¬ 
proximated.  If  the  goal  is  solely  to  perform  simulations,  this  can  be  easily  accomplished  using 
the  same  ODE  solver  used  to  integrate  the  state  equation  (indeed,  the  state  and  estimator 
equations  can  be  combined  into  a  single  system  and  integrated  simultaneously).  This  is  not 
practical  when  experimentally  implementing  the  method,  however,  and  one  must  typically  per¬ 
form  the  work  subject  to  the  following  criteria.  The  method  must  be  sufficiently  efficient  so 
as  to  facilitate  real-time  implementation  and  sufficiently  accurate  so  as  to  resolve  system  dy¬ 
namics.  The  systems  are  quite  often  stiff  which  implies  that  either  a-stability  or  a-stability  is 
important.  Finally,  the  difficulties  in  storing  past  data  make  it  prohibitive  to  use  many  popular 
multistep  methods. 

For  the  experiments  performed  with  the  circular  plate,  the  sample  rate  was  sufficiently  fast 
(and  hence  At  was  sufficiently  small)  that  a  modified  backward  Euler  or  trapezoidal  method 
produced  adequate  results.  We  illustrate  here  such  a  modified  backward  Euler  method. 

4.1.1  Initial  Displacement  and  Velocity 

Considering  first  the  compensator/estimator  sj''stem  with  no  primary  exogenous  force,  we  find 
that  the  modified  Euler  approximate  to  the  solution  at  time  tj^i  is  given  by 

2/cj+i  =  [I  -  AtAo)~^  yoj  +  (/  -  AtAc)~~^  F yob(tj) 

~  A^cUcj  T  F ~yobd’j)  ■ 

The  method  is  modified  in  the  sense  that  current  observation  values  yobiij)  are  used  as  input 
since  futures  values  at  tj+i  are  unknown.  The  time  step  At  is  dictated  by  the  sample  rate.  We 
point  out  that  the  matrix  Ac  =  {I  —  AtAc)~^  and  vector  =  {I  —  AtAc)~^  F  can  be  compnted 
offline  and  then  loaded  as  datafiles  for  the  online  computations.  Hence  the  implicit  nature  of 
the  method,  which  is  necessary  to  ensure  stabilit}^,  does  not  slow  the  implementation.  The 
discrete  time  implementation  of  the  method  is  summarized  in  Algorithm  4.1.1  which  follows. 
The  definitions  of  the  component  matrices  can  be  found  in  Section  3,  and  Algorithm  4.1.1  can 
be  compared  with  the  corresponding  continuous-time  Algorithm  3.1.1  given  in  that  section. 
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Algorithm  4.1.1:  Discrete  Time  Control  of  Initial  Displacement 

and  Velocity 


Offline  (f) 

Construct  matrices  A,  B,C,Q,  i?,  Q,  R 

{ii} 

Solve  Riccati  equations  (3.3)  for  II  and  P 

(Hi) 

Construct  K  =  B^TL 

F  =  PCFR-^ 

A  =  4  -  BK  -  FC 

(iv) 

Choose  appropriate  At  (determined  by  sample  rate) 

iv) 

Construct  Ac  =  (I  —  AtAc)~^ 

Fc  =  (I-  AtAc)~'^  F 

Online  (i) 

Collect  acceleration  data  a{tj) 

(ii) 

Integrate  to  obtain  yoh(tj)  =  v(tj) 

(Hi) 

Time  step  the  discrete  estimator  system 

ycj^-i  Acycj  T  Fcyob(t'j) 

(iv) 

Calculate  the  voltage  u{tj)  =  —Kycj 

4.1.2  Periodic  Exogenous  Force 

The  application  of  these  control  technic|ues  to  systems  with  both  transient  and  steady  state 
behavior  involves  the  approximation  of  both  the  tracking  and  state  estimator  equations  be¬ 
fore  a  control  input  is  calculated.  While  a  variety  of  techniques  and  strategies  can  be  used 
to  obtain  approximate  values  of  r{tj),  which  are  then  used  when  computing  y(tj),  these  cal¬ 
culations  must  ultimately  be  performed  in  real  time  when  implementing  the  method.  In  the 
experiments  involving  the  circular  plate,  the  exogenous  force  was  measured  for  several  periods 
and  the  solutions  to  the  tracking  ecjuation  were  approximated  and  stored  over  a  time  i^eriod 
commensurate  with  the  driving  frecpiency.  These  stored  tracking  values  were  then  used  as  a 
filter  when  approximating  the  estimated  state  during  the  remainder  of  the  experiment. 

Illustrating  with  the  backward  Euler  discretization,  the  approximate  to  the  tracking  solution 
was  determined  from  the  difference  eciuation 

Tj+i  =  (I  —  AtAir)  *  7'j  {I  —  AtAfr)  9i^j) 

=  AiiTj  +  Atrdiij) 

subject  to  the  final  condition  ?’(r)  =  0  (wdien  implementing  the  method,  one  can  simply  search 
for  a  “suitable”  zero  crossing  to  start  the  approximation). 

This  approximation  was  continued  throixghout  several  periods  of  the  driving  force  with  rj 
being  stored  in  a  circular  buffer.  This  buffer  was  then  treated  as  a  filter  when  estimating  the 
state  using  the  difference  equations 

=  (/  -  AtA,)-^  -F  (I  -  AtAc)-^  Fyob(tj)  +  U-  AtAc)-^  A^Vj  F  (/  -  AtA,)-^  git/) 

—  AcVaj  T  F zyobi^j)  T  Ar7 j  -f-  Acg{tj)  . 
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As  indicated  in  the  following  algorithm,  the  time  intensive  calculations  involving  matrix  con¬ 
struction,  solution  of  the  Riccati  equations  and  matrix  inversion  were  performed  offline  prior 
to  the  experiments,  and  the  matrices  Atv)  Ar  and  vector  !Fc  were  simply  loaded  as  datafiles. 
This,  combined  with  the  solution  of  the  tracking  filter  before  state  estimation,  yielded  an  al¬ 
gorithm  which  was  sufficiently  fast  for  implementation.  Current  efforts  are  aimed  toward  the 
simultaneous  approximation  of  the  tracking  and  state  estimator  equations  during  implementa¬ 
tion. 

Algorithm  4.1.2:  Discrete  Time  Control  of  Periodic  Excitation 


Offline 

(*) 

Construct  matrices  A,  B,  C\  Q,R,Q,R 

{ii) 

Solve  Riccati  equations  (3.3)  for  11  and  P 

{in) 

Construct  K  =  R~^B^Il 

F  =  PC^R-^ 

Ac  =  A-  BK  -  FC 

Atr  =  -  [A  -  BKf 

(iv) 

Choose  appropriate  At  (determined  by  sample  rate) 

[v) 

Construct  Ac  =  {I  —  AtAc)~^ 

Ar  =  {I-  AtAc)~^  A, 

Atr  —  {I  ~  At  Air) 

Tc  =  (/  -  AtAc)-^  F 

Online 

(0 

Collect  acceleration  data  a{tj) 

{ii) 

{in) 

Integrate  to  obtain  yoh{tj)  =  v{tj) 

Approximate  and  store  tracking  values 

—  Atj-Vj  1“  Atry{tj) 
r{T)  =  0 

{in) 

Time  step  the  discrete  estimator  system 

2/cj^l  ~  AclJcj  1”  -F ::yob{tj)  1“  AfVj  Acy{tj) 

yc{kT)  =  0 

{iv) 

Calculate  the  voltage  u{tj)  =  —Fycj  +  R~^B^rj 

4.1.3  Higher-Order  Approximations 

In  the  previous  discussion,  a  modified  backward  Euler  method  was  used  to  discretize  the  state 
estimator  and  tracking  ecjuations.  As  indicated,  by  numerical  simulations  reported  in  [14] 
and  experimental  results  in  the  next  section,  for  small  Ai,  this  provides  sufficient  accuracy  to 
calculate  an  effective  feedback  voltage.  If  more  accuracy  is  needed,  a  trapezoid  rule  or  hybrid 
method  of  the  nature  discussed  on  page  225  of  [26]  can  be  used.  These  provide  increased 
accuracy  without  adding  complexity  during  implementation  since  the  components  Atr^Ac^Ar 
and  Tc  can  still  be  computed  offline. 
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4.2  Example  1  -  Matrix  Construction 

The  control  discussion  thus  far  has  been  general  in  the  sense  that  it  holds  for  general  systems 
of  the  form 


ilit)  =  Ay[t)  +  Bu{t)  m  =  Ayit)  +  Buit)  +  g(t) 

or 

yob(t)  =  Chj(t)  Vobit)  =  Cy(t} 

as  long  as  the  pair  {A,  B)  is  stabilizable  and  {A,  C)  is  detectable.  Moreover,  the  cost  functional 
matrices  Q,  R  and  observation  matrices  Q,  R  have  been  treated  as  general  design  criteria  to  be 
specified  according  to  the  application  under  consideration.  In  this  example,  we  illustrate  explic¬ 
itly  the  matrices  and  filters  used  when  implementing  these  control  techniques  for  a  vibrating 
circular  plate. 

We  first  note  that  iV  =  16  modified  cubic  splines  (see  (2.3))  were  sufficient  for  resolving  the 
plate  dynamics  in  the  frec|uency  range  under  consideration.  Due  to  the  axisymmetric  excitation 
and  response  of  the  plate,  the  Fourier  limit  AI  =  0  was  used  in  all  calculations.  Hence  a  total 
of  W”  =  16  basis  functions  were  used  which  led  to  32  coefficients  in  the  vector  y. 

The  formulation  and  sizes  of  all  components  in  the  control  system  for  the  circular  plate  are 
summarized  in  Figure  4.2  and  4.3  below. 


Component 

Size 

Comments 

A  = 

0  /  ' 

32  X  32 

The  elements  comprising  the 
16  X  16  matrices  and 

are  summarized  in  (2.4) 

B  = 

\  i 

'Cq 

o 

1 _ 1 

32  X  1 

See  (2.4)  for  the  description  of 

gH)  = 

\ _ \ 

32  X  1 

The  elements  of  {t)  are  de¬ 

tailed  in  (2.4) 

1  —y 

o  o 

o  o 

<(ri,0i),...,<(rx,0i)  ■ 

p  X  32 

In  the  experiments,  accelera¬ 
tion  data  was  integrated  to 
obtain  velocity  values  which 
are  the  second  state  values  in 
the  second-order  formulation. 
Since  one  accelerometer  was 
used,  p  =  1.  See  Example  2 
of  Section  4.3  for  a  discussion 
regarding  the  duality  between 
control  and  observation. 

Figure  4.2.  System  and  control  matrices  used  in  the  circular  plate  experiments. 


Component 

Size 

Comments 

1 - ( 

o 

to 

o 

1  1 

O 

1 _ 1 

32  X  32 

A  weighted  mass  matrix  was 
used  for  the  penalty  term  Q, 
As  discussed  in  [5],  this  pro¬ 
vides  a  means  of  weighting  the 
kinetic  and  potential  energy  of 
the  plate. 

Example  5.1:  di  =^2  —  1 
Example  5.2:  di—d2  —  b 

R  = 

1 - 1 

Vi 

1 _ 1 

S  X 

In  the  plate  experiments,  one 
controlling  patch  was  used  so 
.s  =  1. 

Example  5.1:  Rn  =  10~^ 
Example  5.2:  Rn  ~  10“^^ 

Q  = 

Ci/^^  0 

0  C2l^ 

32  X  32 

For  the  plate  experiments,  Q 
was  treated  solely  as  a  design 
parameter  as  compared  with 
the  choice  in  [4]  where  physical 
arguments  were  used  to  con¬ 
struct  the  matrix.  The  iden¬ 
tity  weights  were  taken  to  be 
Cl  -=C2  =  1  in  the  experiments. 

R  = 

Rii 

Rpp 

P  X  J) 

In  the  experiments,  p—  \  since 
one  accelerometer  was  used  for 
data  collection.  The  weight 
was  taken  to  be  Rn  =  1. 

Figure  4.3.  Control  and  observation  matrices  used  in  the  cii’cular  plate  experiments.  The  ex- 
irerimental  results  in  Example  5.1  demonstrate  control  of  transient  dynamics  while  the  periodic 
case  is  illustrated  in  Example  5.2. 

These  matrices  and  vectors  are  then  employed  in  Algorithms  4.1.1  or  4.1.2  to  create  the  imple¬ 
mentation  matrices  and  filters  which  were  ultimately  used  in  the  experiments. 
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4.3  Example  2  -  Duality  Between  Control  and  Observation 

The  duality  between  control  and  observation  can  be  noted  by  considering  the  form  of  the  control 
matrix  B  and  observation  matrix  C.  Illustrating  with  the  case  in  which  i  patch  pairs  are  used 
for  control,  the  2M  x  s  control  matrix  has  the  form 

0 

where  the  A'"  x  s  matrix  has  elements 

\b^\  =  / 

L  JTq 

=  f  c,  vW  <h 

J  patch 

(here  denotes  the  characteristic  function  over  the  patch). 

When  data,  from  accelerometers  located  at  the  points  (r^,  ),•••,  (7’p,  is  integrated  to 

obtain  velocitj^  values,  the  p  x  2Af  observation  matrix  is  given  by 


C  =  [0  1] 


c 

0 


where 

hi,.  = 

=  Bf(T„e,). 

With  C  thus  defined,  it  can  immediately  be  noted  that 

yob{t)  =  Cy{i) 

k=i 

denotes  the  physical  value  of  the  velocity  at  the  point  {rj,9j}  given  by  the  state  equations  at 
time  t.  This  is  an  approximation,  to  within  modeling  and  processing  error,  of  the  actual  plate 
plate  velocity  v(t)  which  is  measured  in  experiments. 

Similarly,  multiplication  of  the  state  estimator  coefficients  Pc  by  C  produces  an  c'st  imate  of 
the  velocity  which  is  then  compared  a.gainst  the  measured  plate  velocity  when  into'grating  the 
state  estimator  ecpiation 

Ut)  =  [^  -  y,{t)  +  PC^R-^  [v(t)  -  Cy,it)]  . 

It  should  be  noted  that  as  the  state  estimates  approach  the  measured  plate  values,  the  estimator 
equation  approaches  the  state  equation 

yit)  =  Ayit)  -  BR-^B^Uyc{t) 

which  is  used  to  model  the  j^late  dynamics. 
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4.4  Integration  of  Experimental  Data 

For  the  experiments  involving  the  control  of  circular  plate  vibrations,  data  consisted  of  accel¬ 
eration  measurements  obtained  from  one  or  more  accelerometers  on  the  plate.  It  was  then 
necessary  to  approximately  integrate  this  data  to  obtain  velocity  values  so  as  to  have  a  state 
variable  for  control  calculations.  An  issue  which  turns  out  to  be  crucial  when  approximately 
integrating  experimental  data  concerns  the  robustness  of  the  integrator  with  respect  to  inexact 
initial  conditions  and  DC  gains  or  biases  (added  constants)  in  the  data.  The  inexact  initial  con¬ 
ditions  can  be  due  to  unknown  sj^stem  contributions,  static  shocks  during  system  connections, 
et  cetera.  While  careful  calibration  can  alleviate  some  of  the  uncertainty  in  initial  conditions, 
it  cannot  fully  eliminate  the  problem.  The  problem  of  gains  or  biases  due  to  small  DC  voltages 
in  the  system  can  also  be  minimized  but  never  fully  eliminated.  Hence  an  integrator  which  is 
minimally  affected  by  uncertain  initial  conditions  and  DC  gains  in  the  data  is  crucial  for  success 
when  approximately  integrating  data. 

Here  we  consider  two  techniques  for  approximately  integrating  acceleration  data  to  obtain 
velocity  values  in  accordance  with  the  relation 


v{t)  =  a{t)  . 


Essentially,  the  idea  is  to  replace  the  integration  by  either  the  first-order  differential  eciuation 


V  -f  Du  — 


RC 


(4.1) 


or  the  second-order  ecjuation 

V  Q,v  -(-  (4-2) 

RC 

(see  [21]).  The  design  parameters  ft  and  RC  are  frequency  and  time  constants,  respectively, 
which  are  chosen  so  that  RC  =  1  and  lo  >  6D,  where  uj  is  the  smallest  observed  frecpiency.  For 
solution,  (4.2)  is  written  as  the  first-order  system 


V 

■  -ft 

ft ' 

V 

— 

+ 

RC 

e 

—ft 

0 

e 

0 

z  —  Az  4-  / 


(4.3) 


where  e  =  —  flv.  The  integration  of  (4. .3)  is  subject  to  the  initial  conditions 


■  u(0)  ■ 

’  ^^0  ’ 

.  e(0)  . 

^'0 

(4.4) 


As  will  be  detailed  in  the  subsequent  discussion,  the  first-order  integrator  (4.1)  is  robust 
with  respect  to  inexact  initial  conditions  but  propagates  DC  gains  in  the  data.  In  the  second- 
order  integrator  (4.2),  or  system  (4.3),  both  disturbance  in  initial  conditions  and  DC  gains 
in  acceleration  data,  are  exponentially  attenuated.  Moreover,  the  frequency  response  of  this 
approximator  is  very  close  to  that  in  the  original  signal  for  lo  >  6D.  Hence  this  latter  method 
provides  an  accurate  and  robust  means  of  approximately  integrating  experimental  data. 
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4.4.1  First-Order  Approximate  Integrator 

Here  we  examine  the  proi3erties  of  the  first-order  approximate  integrator  (4.1).  We  consider 
first  the  case  in  which  the  exact  initial  condition  no  is  known  and  the  acceleration  a{t)  is  free 
from  DC  gains  or  biases  (added  constants).  If  we  let  F(s)  =  /i{n(t)}  and  A(s)  =  /i{a(f)},  then 
Laplace  transformation  of  the  sj^stem  (4.1)  yields 


5"  -|-  0 


s  +  n  RC 


Hence  the  approximate  integrator  (4.1)  is  a  single  pole  filter.  Inverse  transformation  then  yields 


v{t)  -  e 


1  f* 


-Q(t-s) 


a{s)ds 


as  the  solution  to  (4.1). 

Similarly,  if  no  denotes  a  perturbed  initial  condition  and  a  DC  gain  %  is  present  in  the  data, 
then  the  solution  is  given  by 


v{t)  =  e  ^*ho 


+  ^  I  ^^[ais)  +  ga]ds 


It  follows  immediately  that 


6(i)  =  ..(i)  +  (0„  -  -  ^) 


Rcn 


It  is  first  noted  that  the  perturbations  in  initial  conditions  exponentially  decay  with  the  rate 
of  decay  influenced  by  the  magnitude  of  the  parameter  0.  DC  gains  of  the  ordeA'  ga/{RCfl) 
remain,  however,  thus  leading  to  difficulties  when  such  biases  are  present  in  the  data.  Both 
properties  are  numerically  illustrated  through  examples  in  [14]. 

The  manner  through  which  the  solution  (4.5)  approximates  the  solution  to  the  original 
relation  v{t)  =  a{t)  can  be  illustrated  with  a  simple  example.  Consider  a{t)  =  TiOTr  cos(1207rt). 
The  solution  to  (4.1)  for  this  acceleration  is 

1  l‘207r  ,  1  120x0  _Q/ 

=  RC  ■  !FT(l2teF  +  120’^  sin(120rt)l  -  + 

which  reduces  to  the  solution  of  the  original  relation  with  0  =  0  and  RC  =  1.  For  0  =  16x, 
the  solution  v  still  provides  an  adequate  approximation  to  the  original  whereas  it  is  a  very  poor 
approximation  with  0  =  120x.  This  phenomenon  is  illustrated  in  [14]. 


4.4.2  Second-Order  Approximate  Integrator 

The  second-order  approximate  integrator  (4.2),  or  equivalent  system  (4.3),  eliminates  the  diffi¬ 
culties  associated  with  both  inexact  initial  conditions  and  DC  gains  or  biases  in  the  data.  The 
elimination  of  constants  in  the  data  can  heuristically  be  attributed  to  the  differentiation  of  the 
acceleration  data.  This  can  be  made  rigorous  by  analytically  solving  the  problem. 
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We  again  let  V{.s)  =  C{v{t)]  and  yl(5)  =  C{a{t)]  and  let  vq.Vi  and  Uq  denote  initial 
conditions.  Transformation  of  (4.2)  yields 

s^V(s)  —  svq  —  Vi  +  fi  [.sy(s)  —  r^o]  +  [5^4(5)  —  aoj 

from  which  it  follows  that 


T//  ^  _  s  +  0,  I  (vi  —  aolRC)  ^ _ s _  1  ^ 

+  ns  +  ^  52  +  ns  +  02  s2  +  Os  +  02  ■  RC  ^  ’  ■ 


Inverse  transformation  then  yields  the  solution 


v{f.) 


cos  ^\/30t/i 


vfflf/2 


^/30t/2 


+  — ^  j  e  cos  ^v^O(t  —  s)/2^  — sin  ^\/30(t  —  .s)/2^  A{s)ds  . 

As  in  the  discussion  of  (4.1),  we  then  consider  the  corresponding  solution  with  perturbed 
initial  conditions  Vo,  Vi  and  So  and  DC  gain  ga-  In  this  case,  the  perturbed  solution  v{t)  is  given 
by 

v{t)  =  v{t)  +  cos  ^\/30t/2^  +  sin  ^■\/30t/2^  [So  —  r>o] 

+  '  tin  [(®‘  ■  m)  -  (”■  -  w). 

,  -at/2 

•/iwn  ■ 

Here  both  the  perturbations  in  initial  conditions  and  the  added  constants  in  the  data  expo¬ 
nentially  decay  with  the  rate  of  decay  dependent  on  the  magnitude  of  D.  This  is  illustrated  in 
examples  given  in  [14].  We  reiterate  that  while  increased  values  of  fl  lead  to  more  rapid  decay  of 
perturbations  and  biases,  the  solution  to  the  differential  equation  less  accurately  approximates 
the  true  velocity.  This,  in  combination  with  the  goal  of  accurately  preserving  signal  frequencies, 
leads  to  the  condition  fl  <  lo/6. 


4.4.3  Numerical  Approximation 


Because  of  the  potential  for  problems  involving  DC  gains  with  the  first-order  formulation  (4.1), 
we  concentrate  primarily  on  the  second-order  filter  (4.2)  and  the  corresponding  system  (4.3). 
We  note  that  similar  scalar  techniques  can  be  used  to  approximate  the  solution  to  (4.1).  In 
considering  numerical  techniques  for  integrating  (4.3),  emphasis  was  placed  on  using  a  technique 
which  could  easily  be  implemented  in  real  time.  The  two  methods  considered  here  are  Euler’s 
method  and  a  backward  Euler’s  method.  The  two  are  summarized  below. 

Euler’s  Method: 

^{tk+i)  —  [^  A  diA]  z{tk)  -1-  dtf(tk) 


1  -  dtn  dtn 
—dtn  1 


diAhl  ' 

RC 

0 


Backward  Euler’s  Method: 


where 


z{tk)  -  [I  -  dtA]  ^  +  [/  -  dtA]  ^  dtf{tk-i) 


[I  -  dtA]-^ 


1 

1  +  dtQ,  +  {dtQ,Y 
—dtil 

1  A  dtO,  A  {dtiiy 


d,tQ 

1  +  dtU  A  (dtfiy 
1  Adt 

1  +  dtAl  A  {dtO.)^ 


dt[I-dtA]  = 


[dtlRC)a[h^^) 

1  +  dtO,  A  (dtfiy 

(-dmiRC)a{h-i) 
1  +  dtfl  A  (dtO,)^ 


The  advantage  of  backward  Euler’s  method  over  Euler’s  method  is  its  stabilit)^  properties 
with  slightly  more  involved  matrices  being  the  disadvantage.  Numerical  examples  demonstrat¬ 
ing  both  methods  with  a  variety  of  exogenous  forces  can  be  found  in  [14]. 


5  Experimental  Results 

Experimental  results  demonstrating  both  the  transient  and  steadj^  state  capabilities  of  the 
control  methodology  are  presented  in  this  section.  The  circular  plate  used  in  these  experiments 
had  a  radius  of  9"  (.2276  ?7?.)  and  a  thickness  of  .05"  (.00127  m).  A  pair  of  piezoceramic  patches 
having  radius  .75"  (.01905  in)  and  thickness  .007"  (.0001778  m)  were  bonded  to  the  center  of  the 
plate.  In  both  the  transient  and  steady  state  experiments,  only  one  patch  was  used  for  control. 
In  the  steady  state  experiments,  the  opposite  patch  was  used  to  drive  the  plate  while  it  was 
allowed  to  remain  uncharged  in  the  transient  case.  The  plate  was  mounted  in  a  wooden  frame 
by  a  circular  aluminum  collar  which  provided  boundary  conditions  which  were  sufficiently  close 
to  clamped  (zero  displacement  and  slope). 

The  first  step  in  the  process  was  the  estimation  of  physical  parameters  through  fit-to-data 
techniques.  As  detailed  in  [2,  3],  transient  plate  vibrations  were  excited  though  an  impact  ham¬ 
mer  strike  or  the  input  of  a  voltage  spike  to  the  patches,  and  acceleration  data  was  measured. 
The  parameter  values  summarized  in  Table  5.1  were  obtained  through  a  least  squares  mini¬ 
mization  of  the  difference  between  the  model  response  and  the  measured  data.  These  values 
were  then  employed  when  constructing  the  component  matrices  A,  B  and  Q  (see  Figure  4.2) 
used  during  the  experimental  implementation  of  the  controller. 
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Physical  Parameters 

/)  •  h 

Plate 

3.170 

{kglnp) 

Plate  +  Pzt 

3.216 

D 

Plate 

11.151 

{N  •  m) 

Plate  +  Pzt 

11.506 

CD 

Plate 

1.443-4 

( N  •  m  •  sec ) 

Plate  +  Pzt 

2.031-4 

u 

Plate 

.326 

Plate  -|-  Pzt 

.325 

CO 

17.021 

IC^(NIV}  - 

Controlling  Patch 

.016 

IC^  (N/V) 

“  Driving  Patch 

.017 

Table  5.1.  Physical  parameters  used  in  the  experiments. 


5.1  Transient  Control 

To  investigate  the  capabilities  of  the  method  for  controlling  transient  vibrations,  decaying  plate 
responses  generated  by  an  impact  hammer  strike  were  considered.  In  each  case,  the  strike  rvas 
directed  to  the  center  of  the  plate  and  hence  the  plate  response  was  axisymmetric.  In  the  first 
set  of  experiments,  data  was  collected  from  an  accelerometer  located  at  the  plate  center  on  the 
side  opposite  from  the  hammer  impact;  thus  P\  =  (ri,^])  =  (0,0)  in  the  construction  of  the 
observation  matrix  C  described  in  Figure  4.2.  The  experiments  were  then  repeated  with  the 
accelerometer  placed  at  the  off-center  point  P2  =  (ri,(9i)  =  (2",  0)  to  illustrate  that  collocation 
between  the  sensor  and  actuator  is  unnecessary  in  this  control  method  (see  Figure  5.1  for 
sensor,  actuator  and  impact  locations).  The  results  obtained  with  the  off-center  accelerometer 
are  presented  here,  and  the  reader  is  referred  to  [14]  for  a  discussion  of  the  transient  control 
results  obtained  with  observations  from  the  centered  accelerometer. 


Figure  5.1.  Patch,  accelerometer  and  impact  locations  for  the  experiments  involving  control 
of  transient  vibrations. 
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Since  no  exogenous  force  was  applied  to  the  plate,  the  state  estimator  and  control  law 
summarized  in  Algorithm  4.1.1  were  used  to  compute  the  controlling  voltage  to  the  patch.  The 
component  matrices  as  well  as  the  cost  functional  and  observation  parameters  used  in  these 
exiseriments  are  summarized  in  Figures  4.2  and  4. -3. 

Data  acquisition  and  processing  were  performed  with  a  PC-based  Texas  Instruments  TMS 
320-C30  digital  signal  processing  (DSP)  board.  A  schematic  of  the  amplifiers,  filters,  DSP 
configuration  and  PC  algorithm  is  given  Figure  5.2.  In  the  experiments,  the  accelerometer 
voltage  was  initially  boosted  by  a  factor  of  10  and  then  reduced  by  5  dB  before  reaching  the 
DSP  board.  The  controlling  voltage  output  from  the  DSP  was  also  boosted  by  an  amplifier 
before  input  to  the  patches.  This  was  necessary  since  the  maximum  voltage  output  by  the  DSP 
is  2.5  V  whereas  60—70  V  were  needed  at  the  patch.  The  reader  is  also  referred  to  Figure  4.1  for 
an  illustration  of  the  experimental  process  and  to  [14]  for  details  regarding  the  implementation 
process.  We  point  out  that  the  ratio  |^,  illustrated  in  the  A/D  conversion,  results  when  the 
voltage  range  -2.5  to  2.5  is  discretized  into  4096  possible  digital  values.  The  reciprocal  process 
occurs  when  digital  values  are  converted  to  analog  voltages  in  the  D/ A  converter. 

The  control  code  was  written  in  assembler  in  order  to  attain  sufficiently  fast  sample  rates 
for  resolving  transient  frequencies  excited  by  the  hammer  impact.  While  the  code  ran  at  rates 
greater  than  7  KHz,  a  sample  rate  of  3.5  KHz  was  used  in  the  experiments.  This  proved  to  be 
sufficient  for  resolving  the  three  axisymmetric  modes  (with  frequencies  of  60  Hz,  227  Hz  and 
512  Hz)  excited  in  the  experiments. 

Representative  plots  of  the  plate  velocity  (integrated  from  the  data)  at  the  off-center  point 
P-i  in  the  uncontrolled  and  controlled  cases  are  given  in  Figure  5.3  with  reduction  levels  at  times 
t  =  0.5, 1.0, 1.5  sec  summarized  in  Table  5.2.  The  percentage  reductions  were  calculated  by  de¬ 
termining  the  ratio  between  the  maximum  values  of  the  controlled  and  uncontrolled  trajectories 
through  one  period  containing  the  time  point  of  interest.  As  illustrated  by  the  results  in  the 
table  and  figure,  the  velocity  level  in  the  controlled  case  has  been  reduced  by  50%  before  .5  sec 
and  is  essentially  fully  attenuated  by  1.5  sec.  We  reiterate  that  these  results  were  obtained  with 
data,  obtained  from  an  accelerometer  at  (2",0)  and  a  centered  actuating  patch  thus  illustrating 
that  collocation  is  unnecessary  for  this  control  method. 

The  voltage,  u(tj)  =  —Kyc^,  was  recorded  in  each  experiment  and  that  yielding  the  control 
results  reported  here  is  plotted  in  Figure  5.4.  It  is  noted  that  the  voltage  has  a  maximum 
magnitude  of  70  V.  In  practice,  it  has  been  observed  that  the  patches  can  be  used  for  extended 
periods  at  the  frequencies  of  interest  without  damage  or  degradation  of  performance  if  the 
voltage  levels  are  maintained  below  8  —  10  rms  V/inil  [28].  Hence  control  voltage  levels  required 
for  control  of  the  transient  vibrations  is  well  within  the  tolerance  of  the  7  mil  patches  used  in 
the  experiments. 

Finally,  the  force  delivered  by  the  hammer  impact  in  the  uncontrolled  and  controlled  cases 
is  plotted  in  Figure  5.5  so  as  to  provide  a  means  of  testing  the  equity  of  excitation  levels  in 
the  uncontrolled  and  controlled  experiments  (the  initial  velocities  in  the  two  cases  can  also  be 
compared  to  determine  whether  the  same  level  of  energy  is  being  delivered  in  each  experiment). 
As  indicated  by  the  results  in  this  latter  figure,  the  force  delivered  in  the  two  cases  is  nearly 
identical  (a  slight  double  hit  was  always  present  when  an  impact  hammer  was  used  to  excite 
the  axisymmetric  modes). 
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Time 

Acceleration 

Velocity 

.5  sec 

69.5 

68.2 

1  sec 

88.9 

84.7 

1.5  sec 

93.8 

97.8 

Table  5.2.  The  percent  reductions  in  acceleration  and  velocity  levels  at  the  point  P2  =  (2",0) 
when  feedback  control  is  implemented. 
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Figure  5.2.  Amplifiers,  DSP  conhguration  and  PC  Algorithm  4.1.1  for  controlling  a  plate 
excited  by  an  initial  impact.  Component  matrices  are  dehnecl  in  Figures  4.2  and  4. -3. 
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Time  (sec) 


Figure  5.3.  Uncontrolled  and  controlled  plate  vibrations  at  (2",  0)  in  response  to  an  impact 
hammer  hit; - (Uncontrolled),  (Controlled). 


Controlling  Voltage  to  the  Patch 


Time  (sec) 


Figure  5.4.  The  controlling  voltage  for  acceleration  data  observed  at  (2",  0). 


Hammer  Forces  in  the  Controlled  and  Uncontrolled  Cases 


Figure  5.5.  The  force  delivered  by  the  impact  hammer  in  the  uncontrolled  and  controlled 
cases, - (Controlled  Case), - (Uncontrolled  Case). 


5.2  Control  of  a  Periodic  Exogenous  Force 

A  second  problem  under  consideration  concerns  the  control  of  plate  vibrations  when  the  plate 
is  driven  by  a  periodic  exogenous  force.  To  demonstrate  the  control  capabilities  in  this  case,  a 
periodic  driving  voltage  was  supplied  to  one  centered  patch  on  the  plate  and  the  patch  on  the 
opposite  side  of  the  plate  was  used  as  the  control  actuator.  Experimental  tests  indicated  that 
a  350  Hz  driving  voltage  produced  a  strong  plate  response  and  all  tests  were  conducted  with 
the  exogenous  voltage  at  that  frequency. 

Two  sources  were  used  to  generate  this  exogenous  signal:  namely,  an  external  oscillator  and 
the  PC  running  the  control  algorithm.  As  reported  under  Case  1  below,  a  purely  steady  state 
response  could  be  considered  with  the  oscillator-generated  exogenous  force  since  the  plate  was 
driven  to  steady  state  before  the  control  program  was  initiated.  Both  a  transient  and  steady 
state  response  were  noted  in  the  PC-generated  signal  since  the  input  of  the  exogenous  force  to 
the  plate  began  at  the  same  time  that  the  control  algorithm  was  started.  This  latter  means  of 
excitation  is  considered  in  Case  2. 

Since  the  system  was  driven  by  a  periodic  exogenous  force,  the  discrete  time  Algorithm 

4.1.2  (corresponding  to  continuous  time  Algorithm  3.3.1)  was  used  to  calcula-te  the  controlling 
voltage  to  the  actuating  patch.  Again,  component  matrices  as  well  as  control  and  observation 
parameters  are  summarized  in  Figures  4.2  and  4.3. 

For  the  off-center  results  reported  here,  velocity  observation  values  were  obtained  by  inte¬ 
grating  data  from  an  off-center  accelerometer  located  at  the  point  P2  =  (2",  0)  as  depicted  in 
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Figure  5.1.  The  second  input  to  the  algorithm  consisted  of  measurements  g{tj)  of  the  force 
driving  the  plate.  We  point  out  that  when  implementing  this  control  method,  both  phase 
and  magnitude  information  for  the  driving  force  were  required  as  compared  with  many  other 
methods  (e.g.,  feedforward)  which  require  only  phase  information. 

In  these  experiments,  the  tracking  components  Vj  were  calculated  first  and  stored  in  a  cir¬ 
cular  buffer.  These  values  were  then  used  when  calculating  the  state  estimates  ijj  and  voltages 
u{tj).  While  the  implementation  in  this  manner  lacilitated  running  the  algorithm  with  sample 
rates  on  the  order  of  7  KHz  (again,  the  algorithm  was  coded  in  assembler),  it  limits  the  robust¬ 
ness  of  the  method  with  respect  to  changes  and  variations  in  the  di’iving  force.  Current  efforts 
are  directed  toward  simultaneous  solution  of  the  tracking  and  estimator  difference  equations. 

A  crucial  issue  when  implementing  the  method  concerns  the  handling  ol  delays  and  phase 
shifts  produced  by  the  filters,  A/D  and  D/A  conversions,  computation  of  the  control  voltages, 
et  cetera.  While  the  amount  of  delay  and  phase  shift  is  frequency  dependent,  experiments 
indicated  that  at  350  Hz,  30  -  40^  phase  shifts  were  introduced  by  the  hardware.  This  was  suT 
ficient  to  destabilize  the  controller  if  left  uncompensated.  In  the  experiments,  we  compensated 
by  first  conducting  an  offline,  numerical  identification’  to  determine  the  amount  of  added  delay 
necessary  for  stabilizing  the  controller  in  the  presence  of  phase  shifts  of  the  order  introduced  by 
the  experimental  hardware.  A  summary  of  these  results  can  be  found  in  [14].  This  is  analogous 
to  the  online  tuning  or  phase  locking  which  is  necessary  for  ensuring  stability  in  other  control 
methods.  The  numerical  tests,  summarized  in  [14],  indicated  that  the  introduction  of  a  216*" 
delay  in  accelerometer  or  exogenous  force  data  would  stabilize  the  systems,  and  this  was  im¬ 
plemented  in  the  experiments  b)^  using  a  DAC  delay  box.  The  compensation  for  phase  delays 
in  this  manner  provides  merely  a  first  step  toward  optimal  implementation  of  the  method  and 
one  aspect  of  current  research  efforts  is  directed  toward  online  compensation  in  the  algorithm. 

Data  accjuisition  and  processing  was  again  performed  with  a  PC-based  Texas  Instruments 
TMS  320-C30  board.  A  schematic  of  the  setup  is  given  in  Figure  5.6. 
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Figure  5. 6-  Amplifiers,  DSP  configuration  and  PC  Algorithm  4.1.2  for  controlling  a  plate 
driven  by  a  periodic  exogenous  force.  Component  matrices  are  defined  in  Figures  4.2  and  4.3. 


Case  1:  Oscillator-Generated  Driving  Signal 

For  the  results  described  here,  the  signal  to  the  driving  patch  was  generated  by  an  external 
oscillator.  The  plate  was  allowed  to  reach  steadj^  state  and  then  the  control  program  was  initi¬ 
ated.  Acceleration  levels  measured  by  the  accelerometer  located  at  P2  =  (2",0)  and  integrated 
velocity  values  for  the  uncontrolled  and  controlled  cases  are  plotted  in  Figure  5.7.  As  noted 
from  the  controlled  trajectories  in  that  figure,  it  takes  the  algorithm  approximately  0.06  seconds 
to  calculate  and  store  a  sufficient  number  of  tracking  components  Vj.  During  that  time  interval, 
no  voltage  is  fed  to  the  actuating  patch.  Once  the  tracking  calculations  are  completed,  state 
estimation  begins  and  the  controlling  voltage  is  computed  and  fed  back  into  the  system.  The 
vibration  levels  decay  for  approximateh^  0.3  seconds  and  then  are  maintained  at  levels  that 
are  approximately  15%  of  those  for  the  for  the  uncontrolled  case  for  the  remaiink'r  of  the  time 
interval.  This  corresponds  to  a  20log{acon/ <^uncon)  ~  —16.5  dB  reduction  in  acceleration  levels. 

While  the  magnitude  of  the  controlling  voltage  is  dependent  upon  the  amplitude  of  the 
driving  signal,  magnitudes  less  than  40  V  (28.3  V  rms)  were  required  to  attain  the  levels  reported 
here.  At  350  Hz,  this  was  well  within  the  range  (56  —  70  V  rms  in  this  case)  that  was  considered 
to  be  safe  for  the  patch  being  used  (see  Example  5.1  for  further  discussion  regarding  the  voltage 
levels  at  which  the  patches  can  be  driven  without  damage  or  degradation  in  performance). 

The  results  in  this  experiment  demonstrate  the  effectiveness  of  the  algorithm  for  controlling 
a  system  that  has  reached  steady  state.  Hence  one  is  not  required  to  start  the  control  with  a 
system  at  rest. 
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Figure  5.7.  Uncontrolled  and  controlled  plate  acceleration  and  velocity  at  (2,0)  for  system 
driven  by  periodic  exogenous  force; - (Uncontrolled),  —  (Controlled). 

Case  2:  PC-Generated  Driving  Signal 

A  second  mechanism  for  generating  the  driving  signal  is  with  the  PC  that  is  used  to  process 
data  and  run  the  control  algorithm.  Acceleration  and  velocity  plots  of  the  uncontrolled  and 
controlled  plate  vibrations  excited  in  this  manner  are  given  in  Figure  5.8.  It  can  be  seen  that  in 
this  case,  the  plate  starts  from  rest  and  is  still  being  driven  through  a  transient  stage  when  the 
tracking  calculations  are  completed  and  control  begins.  At  that  point,  the  controlled  trajecto¬ 
ries  are  reduced  to  the  levels  noted  in  the  purely  steady  state  case  whereas  the  uncontrolled 
trajectories  are  driven  to  steady  state.  Here,  an  82%  (15  dB)  reduction  in  levels  is  noted  at  time 
T  =  I  second.  This  was  obtained  with  a  controlling  voltage  of  magnitude  12  (8.5  U, -ms  ). 

These  results  demonstrate  the  effectiveness  of  the  control  algorithm  for  a  system  undergoing 
transient  oscillations  before  reaching  steady  state  in  response  to  a  periodic  driving  force. 


Uncontrolled  and  Controlled  Plate  Acceleration  x  10'®  Uncontrolled  and  Controlled  Plate  Velocity 


Time  (sec)  Time  (sec) 

Figure  5.8.  Uncontrolled  and  controlled  plate  acceleration  and  velocity  at  (2,0)  for  system 
driven  by  periodic  exogenous  force; - (Uncontrolled),  (Controlled). 


6  Conclusions 


In  this  work,  the  experimental  implementation  of  a  PDE-based  controller  was  considered.  While 
the  motivating  application  involves  the  control  of  vibration  levels  for  a  circular  plate  through 
the  excitation  of  surface-mounted  piezoceramic  patches,  the  general  technicpies  desci'ibed  here 
will  extend  to  a  variety  of  applications. 

For  such  control  techniques,  the  first  step  is  the  derivation  of  a  PDE  model  which  accurately 
describes  the  dynamics  of  the  system  under  consideration.  Physical  parameters  in  these  models 
must  typically  be  estimated  through  fit-to-data  techniques  before  control  applications  can  be 
considered.  Following  a  brief  discussion  regarding  the  strong  and  weak  forms  of  a  thin  plate 
model  with  a  discontinuous  control  input  term  (due  to  the  piecewise  constant  nature  of  the 
piezoceramic  patches),  continuous-time  LQG  and  methods  for  systems  with  no  exogenous 
force  or  a  periodic  exogenous  force  were  discussed.  The  discrete-time  approximations  necessary 
for  implementing  the  methods  with  digital  measurements  were  also  pre,sented.  A  crucial  step 
when  implementing  the  discrete-time  controllers  involves  the  approximate  integration  of  data 
(e.g.,  accelerometer  data  integrated  to  obtain  velocity  state  values),  and  first  and  second-order 
filters  for  accomplishing  this  were  discussed.  Without  such  filters,  DC  biases,  that  are  always 
present  in  the  data,  would  render  the  integrated  values  useless. 

Experimental  results  demonstrating  the  control  of  transient  and  steady  state  vibrations  were 
then  presented.  One  advantage  of  the  PDE-based  controllers  over  standard  frecpiency  response 
input/output  techniques  is  the  capability  for  direct  control  of  transient  responses,  and  this  was 
demonstrated  in  the  first  set  of  examples.  A  centered  hammer  impact  was  used  to  excite  the 
plate  and  integrated  data  from  an  off-center  accelerometer  was  used  to  reconstruct  the  state. 
The  results  demonstrate  that  attenuation  levels  on  the  order  of  95%  reduction  can  be  attained 
by  1.5  seconds  using  the  PDE-based  controller. 

The  second  example  demonstrates  the  control  of  transient  and  steadj^  state  responses  when 
the  plate  was  driven  by  a  periodic  exogenous  voltage  to  a  secondary  piezoceramic  patch.  These 
results  demonstrate  that,  after  accounting  for  hardware  delays,  attenuation  levels  on  the  order 
of  85%  were  attained  when  control  was  implemented.  While  implementation  techniques  are 
still  being  refined,  these  results  demonstrate  the  effectiveness  of  the  PDE-based  controller  for 
this  system  and  indicate  the  potential  of  these  control  techniques  for  reducing  transient  and 
steady  state  dynamics  in  other  structural  and  structural  acoustic  sj^stems. 
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